Bismuth in strong magnetic fields: unconventional Zeeman coupling and correlation effects 



00 

o 
o 

(N 
•*-> 

o 
O 
o 



U 

a 

I 

c 

o 
o 



(N 

m 
d 

oo 
o 



Jason Alicea 1 and Leon Balents 2 

'Department of Physics, California Institute of Technology, Pasadena, CA 91125 
2 Kavli Insitute for Theoretical Physics, University of California, Santa Barbara, CA 93106-9530 

(Dated: October 20, 2008) 

Recent experiments on bismuth have uncovered remarkably rich magnetization structure at fields well beyond 
the regime in which all carriers are expected to reside in the lowest Landau level. Motivated by these findings, we 
start from a microscopic tight-binding model and derive a low-energy Hamiltonian for the holes and three Dirac 
electrons pockets in bismuth. We find that an unconventional electron Zeeman effect, overlooked previously, 
suppresses the quantum limit for the electrons dramatically, giving rise to the observed anomalous magnetization 
structure. We further study interaction effects near fields at which the 2nd Landau level for one electron pocket 
empties, where magnetization hysteresis was observed. Here we find instabilities towards both charge density 
wave and Wigner crystal phases, and propose that hysteresis arises from a first-order transition out of the latter. 

PACS numbers: 



Introduction. After more than a century of active research, 
bismuth continues to yield fascinating discoveries. Much of 
this material's exceptional behavior stems from its band struc- 
ture HI — the Fermi surface arises from a hole pocket and three 
Dirac electron pockets which contribute an extremely low car- 
rier density of 3 x 10 18 cm~ 3 . One remarkable consequence of 
the small carrier concentration is that weak Sb doping is be- 
lieved to change the system from a semimetal to a topological 
insulator yl |3fl. Another is that the quantum limit, at which 
carriers are confined to the lowest Landau level (LLL), can be 
realized with laboratory fields. In particular, for fields along 
the highest-symmetry 'trigonal' axis, the hole quantum limit 
occurs at ^9T |4j[. The electron quantum limit, while less 
clear experimentally, is believed to occur at similar fields fl. 

From a single-particle perspective, quantum oscillations 
should subside once all carriers reside in the LLL, and trans- 
port and thermodynamic quantities should appear 'feature- 
less'. Recent experiments on bismuth in trigonal fields nev- 
ertheless observed surprisingly rich physics extending well 
beyond 9T. Nernst peaks were resolved at 13.3, 22.3, and 
30. 8T, with the Hall resistance exhibiting step-like features in 
between, prompting the suggestion that fractionalization may 
be occurring in this 3D system [6]. More recent torque mag- 
netometry studies at fields near the trigonal axis additionally 
measured unanticipated magnetization structure — including 
hysteresis — persisting up to the largest fields probed (3 IT), 
providing further evidence for correlation physics Q3Q . 

Here we derive a low-energy theory for bismuth, suitable 
for analyzing the Landau level (LL) structure and interaction 
effects. For the electrons, we show that strong spin-orbit cou- 
pling generates unconventional Zeeman terms not present in 
the standard Cohen-Blount model [7]. This coupling sup- 
presses the electron quantum limit far beyond 9T for trigo- 
nal fields contrary to what has been assumed, recovering the 
high-field magnetization structure observed in Ref. 5 [see Fig. 
[Ha)]. Turning to Coulomb effects, we show that interactions 
have the strongest influence near fields for which a low LL 
empties. In agreement with this expectation, the experimen- 
tally observed hysteresis coincides with one electron pocket 
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FIG. 1: (a) Single-particle DOS, excluding the contribution from 
the electron pocket invisible in torque experiments All electron 
pockets occupy the 2nd LL in region I, while pocket 3 empties into 
the LLL in region II. (b) Schematic energy dispersion for pocket 3 
and (c) proposed phase diagram with interactions near B*(9). 



emptying its 2nd LL. To address this aspect of experiment, we 
study interaction effects near this band emptying, employing 
Functional Renormalization Group (FRG) techniques similar 
to Ref. 8 and going beyond the early analysis of Abrikosov 
Q9[] . As the Fermi energy passes from above the 2nd LL into 
the LLL, we find that the leading instability for these electrons 
involves charge-density- wave (CDW) order in the 2nd LL, fol- 
lowed by Wigner crystal formation and LLL CDW order [see 
Fig. [He)]. We suggest that the hysteresis originates from the 
Wigner crystal phase, and discuss experiments to verify this 
proposal. Finally, we comment briefly on the Nernst and Hall 
effect puzzles, which remain unexplained by this work. 

Low-energy theory. Our starting point is the tight-binding 
model of Liu and Allen ifioll . which was constructed to accu- 
rately reproduce bismuth's band structure near the Fermi level 
Jjl. We derive from this microscopic model an effective low- 
energy Hamiltonian for the electrons and holes in a magnetic 
field B by expanding the lattice fermion operators as follows, 
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where fi — s,p XtV . z labels the outer-shell orbital with spin 
a — T / I at s i te r on sublattice j = 1,2. The two-component 
operators h describe hole excitations near wavevector Q and 
the four-component Dirac fermions ip\ describe electron exci- 
tations near P>; the corresponding wavefunctions are ip° and 
<p xe . Using Eq. ([T), we obtain the non-interacting Hamiltonian 
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Here D = V iqA, with B = V x A and q = ±e 
for the holes and electrons, respectively. We employ three 



sets of Pauli matrices a 

rf = fi x (v lx T y 
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nians Hi' 2 for pockets 1 and 2 can be obtained from H\ 
by rotating D,B by ±27r/3 about the trigonal (z) axis. In 
terms of the electron mass m e and speed of light c, we 
have fi e = 0.0335eV, fx h = 0.012eV, m ± = 0.0675m e , 
m z = 0.612m e , m = 6meV, u la! = 0.0022c, v 2x = -0.002c, 
uij, = -7.6 x 10~ 5 c, w 2y = 3.4 x 10" 4 c, v lz = 0.002c, 
V2z — —0.0014c. The Fermi energy Ep — when B = 0, 
but changes to maintain charge neutrality when B^O. 

Zeeman coupling warrants further discussion. Microscop- 
ically, Zeeman energy has one source oc B • (L + 2S) and 
another coming from spin-orbit coupling oc (VU x A) • S, 
where U is the crystal potential and L, S denote orbital/spin 
angular momentum. While the corresponding low-energy 
terms can be obtained via Eq. ([T), evaluating the spin-orbit 
contribution is nontrivial since the potential U is unknown. 
Progress can be made, however, by assuming U is rotation- 
ally invariant, as then only two independent matrix elements 
remain: x s / P = T^^{ s /Px T \rd r U sin 2 9\s/p x f), where 
\s/p x |) are the atomic s/p x orbital wavefunctions. With 
this assumption, we obtain the hole Zeeman splitting in Eq. 
(O, which is sensitive only to B z , consistent with exper- 
iment [4J. For the electrons, we obtain the Zeeman cou- 
pling in Eq. © with G x = r x (g lx + g 2x ^ z ) and G VtZ = 
T v {gi v ,z + g2y,zH z ) + r z (g 3y , z + 54y,zM z )- The g-factors 
are listed in Table U Crucially, gh depends only on Xp, and 
the electron g-factors obtain only a weak contribution from 
Xs because the wavefunctions <ft xe are concentrated on the p- 
orbitals. The precise value of Xs is therefore unimportant, 
and we will simply set Xs = Xp- Finally, hole Zeeman split- 
ting has been well-studied experimentally 111 111 , and from the 
available data we deduce that gh 54, which allows us to 
determine Xv anc ^ nence tne electron q-factors. lfl7ll 

Equations (fJI and (f3]l constitute the first major result of this 
paper. Most importantly, the electron Zeeman coupling has 
not been considered previously, and modifies the spectrum 
dramatically at high fields as we now discuss. 

Single-particle spectrum. While the hole Hamiltonian is 
easily diagonalized, the electron part is much more difficult 
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TABLE I: Electron/hole g-factors, with g a — a a + b a \v + c «Xs 



since the Zeeman and orbital terms do not commute. Qual- 
itatively, Zeeman coupling generates components of higher 
LLs in the wavefunctions compared to the orbital-only prob- 
lem and splits the usual LL degeneracy in the Dirac spectrum 
(see, e.g., Supporting Online Material for Ref. 5). To proceed, 
we truncate the Hilbert space, including only the first n <~ 10 
LLs, and diagonalize the Hamiltonian numerically. 

Motivated by magnetization experiments [5], we have stud- 
ied the spectrum for fields tilted by an angle from the trig- 
onal towards the binary axis. Figure [T(a)| displays the single- 
particle density of states (DOS) in the B — 9 plane (excluding 
the electron pocket invisible to experiments (5[]). Bright lines 
occur where the Fermi energy crosses the bottom of a LL, and 
the flat line at ~9T corresponds to the hole quantum limit. 
Remarkably, the electrons give rise to features at much higher 
fields which agree well with the anomalous structure reported 
experimentally (see Fig. 3 in Ref. [5b. 

The persistence of this structure up to such large fields 
arises from an increase in carrier density with field ||4j] and 
the splitting of the electron LL degeneracy by Zeeman cou- 
pling, which makes 2nd LL states available at lower energies 
and strongly suppresses the electron quantum limit. Indeed, in 
region I of Fig. 1 1 (a)| all electron pockets occupy the 2nd LL. 
Tilting the field into region II pushes one of those pockets (as- 
sumed to be pocket 3 hereafter) into the LLL. The dispersion 
for pocket 3 versus the momentum fcii along the field appears 
schematically in Fig. [TJb), together with the chemical poten- 
tial in I and II. Interestingly, experiments observed hysteresis 
in the magnetization at the line labeled (9) separating these 
regions [corresponding to //* in Fig. [TJb)]. Addressing this 
puzzle requires moving beyond single-particle physics. 

Correlation effects. We now add Coulomb interactions, 
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where p = J^x tyypx — f^h and V(r) is the screened 
Coulomb potential. Note that the number of holes and elec- 
trons within each pocket remains separately conserved here. 
Short-range pieces arising from lattice effects break this sym- 
metry, but are subdominant and can be neglected. 

Interactions can be most simply treated at weak coupling, 
which is controlled provided the dimensionless interaction 
strengths satisfy e 2 /w™ <C 1, where u" is the Fermi velocity 
for pocket a in LL n. This criterion inevitably breaks down 
near a band emptying since at least one u™ — > 0. Correla- 
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tion effects will be most pronounced in this strongly coupled 
regime, and this is indeed where hysteresis appears. Our aim 
now is to study the crossover from weak to strong coupling to 
understand this aspect of experiment. Rather than addressing 
the fate of all carriers, we will focus more narrowly on the 
leading instability involving pocket 3 electrons near B*(9). 

Even this restricted question requires considering numer- 
ous interactions, since H mt couples pocket 3 to all other car- 
riers, some of which occupy more than one LL. We can, how- 
ever, further distill the problem using an intuitive principle: 
similar carriers generally couple more effectively than dis- 
similar ones. For instance, velocity anisotropies sharply dis- 
tort the electron LL wavefunctions, which strongly suppresses 
instabilities involving different electron pockets. [18] Further- 
more, the hole and electron Fermi velocities differ signifi- 
cantly, which disfavors electron-hole pairing. We will there- 
fore concentrate on interactions involving only pocket 3 elec- 
trons, commenting briefly on other couplings. For the three 
important cases [regions I, II, and the line B* (0)], we employ 
the FRG method to determine the leading instability. FRG 
equations give the renormalized interactions at length scale L 
as a function of the logarithmic rescaling factor I = ln(L/A), 
where A is a microscopic length of order the Fermi wave- 
length. 

We begin in region II, where pocket 3 is confined to the 
LLL. To obtain a low-energy Hamiltonian we linearize the ki- 
netic energy around the Fermi momenta ±fc^ 3 , expand ^3 in 
terms of right/left movers, and project onto the LLL (we em- 
ploy Landau gauge and label the transverse momentum k±). 
Interactions between pocket 3 electrons can then be written 

^int" = J p{k±i,k_ L2 )c°^ 3 (k 1 +k 3 )c° t 3 (k 2 +k 3 ) 

x C° 3 (k 3 )4 3 (k 1+ k 2 + k 3 ), (5) 

where k = (k±,k^) and c°^ L3 creates a right/left-moving 
LLL electron. The FRG equation for the coupling function 
p(k, x) = p(r) Fourier transformed in the second argument is 
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with (f*g)(r) = f r , r „e^ rAr ' +r ' Ar " +r " / ^f(r')g(r"). The 
solution to Eq. (O is well-understood H8L Hal . The first term 
in the (3 function causes p to flow off, driving condensation 
of (c^ R c® L ) which gaps this channel and yields 2k° F3i CDW 
order along the field; the second term merely reduces the tran- 
sition temperature [8]. Furthermore, it follows from the initial 
conditions that p(r = 0) is largest at the instability, so that the 
transverse density is uniform in the CDW phase. As an aside, 
we note that coupling to other electron pockets and the holes 
leads to a series of couplings with flow equations analogous to 
Eq. ©. Due to the velocity anisotropy discussed above, inter- 
pocket electron instabilities set in at very low energies, and are 
preempted by the p instability. Similarly, since the hole veloc- 
ity is several times smaller than v® , hole-hole pairing preempts 
the 'exci tonic' instability [13] in the electron-hole channel. 



Next, we sit at (0) and fine-tune the chemical potential to 
/i» in Fig. 01b), precisely at the bottom of the 2nd LL. At low 
energies, we now have linearly-dispersing right/left movers 
C °R,/L3 cou pl e d w i m quadratically dispersing 2nd LL electrons 
which we denote by d\. As it stands, the problem can not be 
treated within weak coupling since the soft dispersion for the 
latter renders interactions involving d\ strongly relevant at the 
non-interacting fixed point. To proceed, we follow Ref. 
and seek a controlled limit by replacing the 2nd LL kinetic 
energy Dk 2 with D|fc|| | 1+e and performing an expansion in e 
and the interaction strength in the limit 1 ^ e ~ e 2 /u 3 . 

Pocket 3 interactions now involve several couplings: 



f • = J {P(k±i,k ±2 )c° 3 \ i (k 1 + k 3 )c%(k 2 + k 3 ) 

4 3 (k 3 )4 3 (k 1 + k 2 + k 3 ) + ud\U\U\d\ <7s 

Pt ot ,1,1 



+ Kk4 T 44 3 + (R - L)} + [wc^cZd'A + h.c.}}. 

The arguments in the last three couplings have been sup- 
pressed for brevity, but should appear as in the p term. Fourier 
transforming in the second argument and defining r = (fc, x) 
as before, it will be convenient below to write the flow equa- 
tions for these coupling functions as follows: 

9ep(r) = ^^[p(r) 2 -(p*p){r)]-?-(w*w*)(r), 

(3 a 
deu(r) — eu(r) H — u(r) 2 (u*u)(r) 

TT TT 

- (8) 
d e v(r) = 2L[ v (rf ~(v*v)(r)] + ^\w(r)\ 2 , 

Z7T 7T 

el ex 
d e w{r) = -w(r) - -^-^(p*w)(r) - -(w*u)(r) 

+ ^-v(r)w(r) - - [ e irAr 'v(r')w(r'). 

TT TT J r , 

Here a = 1/D,7i = $= l/{vl+DX e e- e % 72 = 7iA £ e" rf , 
and (3 = 0. Naively, power-counting suggests that u flows off 
first. However, since (3 = 0, the analogous term that drove 
the CDW instability in region II is absent, which leaves open 
a more interesting possibility. 

To facilitate analytic and numeric progress, we now approx- 
imate V(r) as local (anisotropies can then be scaled away). 
Note first that at the initial conditions, w is suppressed com- 
pared to p, u, v since terms that survive Pauli exclusion are 
suppressed due to small overlaps between the participating LL 
wavefunctions. Thus we begin by setting w = 0, which de- 
couples the remaining equations. Clearly both p and v then 
flow off, with the former diverging faster since p is larger at 
the initial conditions. The behavior of u is less obvious. As- 
suming u(k,x) is rotationally invariant (after rescaling) and 
analytic in k,x, we can solve the flows for u by assuming an 
ansatz it(r; t) = J2j fj(^)Xj( r )- Here j runs over 0, 2, 4, . . . 
and Xj( r ) — e ~ r2 ^ 2 Pj( r )> with Pj{ r ) degree-j polynomials 
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defined so that (xj * Xj')( r ) — dj,j'CjXj( r )- One then ob- 
tains decoupled equations for each fj which show that u(r; £) 
flows to a finite fixed point provided Cjfj (0) > for all j. In 
our problem the initial conditions are well-approximated by 
u(r;£ — 0) = c(r 2 — l)e~ r I 2 cx Xiif), with c > 0, which 
indeed remains finite under renormalization. Though this re- 
sult is correct in the e-expansion, we caution that the fixed 
point could move to strong coupling for physical e = 1. This 
would imply a it-driven instability (see below). 

To treat the full problem with w ^ we rely on numer- 
ical integration of the FRG equations. Here we find that w 
further enhances the instability in p, and also generates addi- 
tional non-zero components of Xj m u , causing this function 
to flow off in tandem with w. In contrast, v is only weakly 
affected by w and remains unimportant. Because p was al- 
ready unstable when w = 0, this function becomes of order 
one first, driving CDW order in the LLL with uniform trans- 
verse density. The leading instability in the 2nd LL occurs 
soon after driven by u, which localizes these states produc- 
ing Wigner crystallization as the chemical potential increases. 
While this is certainly reasonable, we caution that definitively 
ruling out the scenario where u provides the leading instability 
is difficult due numerical limitations, particularly when V(r) 
is non-local. In this case the order of the instabilities is simply 
reversed. 

Finally, consider region I, with the chemical potential 
slightly above the 2nd LL for pocket 3. This problem is rem- 
iniscent of a ID wire with 2 transverse modes occupied 1151 . 
though the LLs change the physics qualitatively. As in region 
II, we derive a low-energy theory for right/left movers in the 
LLL, c^ R , L , and 2nd LL, c\ R , L . Interactions then read 

H?£ = J {p(fcxi, fcl2)c&(k! + k 3 )c%(k 2 + k 3 ) 

xc°3(k 3 )4 3 (k 1 + k 2 + k 3 ) + uc^ R c\lcl L cl R (9) 

+ [ vc °3R c Il c Il c R3 + (R^ L)] + [wcf R c^ L c\ L c\ R + h.c.]}, 

where again the arguments of the last three couplings should 
appear as in the first. These couplings flow as in Eqs. (0, 
but with e = 5 = 0, a = /3 = l/(2wg), and 71 = 472 = 
2/(wg + U3). It is natural to suspect here that u dominates, 
since the 2nd LL carriers have the slowest velocity. This is 
indeed correct, which can be understood analytically by ig- 
noring all (/ * g) terms in Eqs. ([8]). Since (3 is the largest 
coefficient remaining when the 2nd LL is weakly populated, 
u flows off before all other couplings, driving CDW order 
formed by condensing {c^„c\ L ). The transverse density be- 
comes Wigner crystalline [ 12] at the instability, since the form 
of the 2nd LL Dirac wavefunctions dictates that ti(r) is max- 
imized at r ^ 0. In the limit where V(r) is local, we have 
verified numerically that terms neglected in this crude analy- 
sis do not modify these conclusions. 

Discussion. Putting our results together, we propose the 
minimal phase diagram shown in Fig. [He). We speculate that 
experiments of Ref. 5 may be conducted above the critical 
temperature for the CDW phases in regions I and II but below 



that for the Wigner crystal phase near the boundary. The lat- 
ter transitions are almost certainly first order, which would be 
consistent with the observed hysteresis. However, two transi- 
tions ought to occur here, whereas experiments see only one. 
This issue can be resolved by invoking disorder, which will be 
particularly important in the low-density region of the Wigner 
crystal phase close to region II. The transition on that side 
is expected to be smeared by disorder-pinning of the local- 
ized states, which should be addressed in more detail. Future 
experiments, particularly nonlinear transport and x-ray scat- 
tering studies to search for signs of CDW and Wigner crystal 
order, should provide valuable clues as to the true nature of 
this transition. 

The Hall and Nernst effect puzzles in trigonal fields B > 
9T ||6|] are difficult to resolve at weak coupling. Here hole-hole 
pairing drives the leading CDW instability, which is not ex- 
pected to recover these anomalies. While surface states should 
be seriously considered as emphasized in Ref. 5, we believe a 
more exotic origin is not inconceivable. Interactions between 
the LLL holes are not weak, as e 2 /v® ~ 24, well outside of 
the range where weak coupling is expected to be reliable. In 
contrast, interactions between LLL electrons, while not nec- 
essarily weak, are several times smaller: e 2 /v® ~ 5. Al- 
though screening should reduce these somewhat, the problem 
warrants studying from a strong-coupling standpoint, which 
presents an exciting research direction 1 16]. 

It is a pleasure to acknowledge illuminating discussions 
with K. Behnia, Lu Li, O. Motrunich, N.R Ong, and G. Re- 
fael, as well as the hospitality of the KITP where this work 
was initiated. We also acknowledge support from the Lee A. 
DuBridge Foundation (JA), and the Packard Foundation and 
National Science Foundation through grants DMR-0804564 
and PHY05-5 1164 (LB). 



[1] V. S. Edel'man, Adv. Phys. 25, 555 (1976). 

[2] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

[3] D. Hsieh et al., Nature 452, 970 (2008). 

[4] G. E. Smith et al, Phys. Rev. 135, Al 118 (1964). 

[5] L. Li et al., Science 321, 547 (2008). 

[6] K. Behnia et al, Science 317, 1729 (2007). 

[7] M. H. Cohen and E. I. Blount, Phil. Mag. 5, 1 15 (1960). 

[8] V. M. Yakovenko, Phys. Rev. B 47, 8851 (1993). 

[9] A. A. Abrikosov, J. Low Temp. Phys. 10, 3 (1973). 
[10] Y. Liu and R. E. Allen, Phys. Rev. B 52, 1566 (1995). 
[1 1] S. G. Bompadre et al, Phys. Rev. B 64, 073103 (2001). 
[12] S.-W. Tsai et al, Phys. Rev. B 65, 241 102 (2002). 
[13] B. I. Halperin and T. M. Rice, Rev. Mod. Phys. 40, 755 (1968). 
[14] L. Balents and M. P. A. Fisher, Phys. Rev. B 53, 12133 (1996). 
[15] J. S. Meyer et al, Phys. Rev. Lett. 98, 126404 (2007). 
[16] F. J. Burnell et al, arXiv:0810.1757 (unpublished). 
[17] The data of Ref. Qj] appear consistent with the hole Zeeman 
splitting being 2 ± 0.16 times the cyclotron energy. We use the 
lower sign, which gives gu = 54, since the hole quantum limit 
then occurs closer to 9T as observed experimentally. 
[18] Fields oriented so that 2 or 3 pockets are symmetry-related may 
produce instabilities involving multiple electron pockets. 



